Do low surface brightness galaxies have dense disks? 
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The disk masses of four low surface brightness galaxies (LSB) were estimated 
using marginal gravitational stability criterion and the stellar velocity dispersion 
data which were taken from Pizzella et al., 2008 111. The constructed mass models 
appear to be close to the models of maximal disk. The results show that the disks 
of LSB galaxies may be significantly more massive than it is usually accepted from 
their brightnesses. In this case their surface densities and masses appear to be rather 
typical for normal spirals. Otherwise, unlike the disks of many spiral galaxies, the 

q LSB disks are dynamically overheated. 

U 

O I- INTRODUCTION 

LSB galaxies represent the wide class of disk galaxies, central surface brightnesses of 

y—i which are many times lower than those of normal spiral and SO galaxies. Among LSB 

> 

galaxies there are both dwarf and giant systems, both objects with solar metallicity and 
systems with significant deficiency of metals. The gas mass over the luminosity Mhi/Lb 
\q ratio of LSB galaxies lies in the wide interval: 0,1-50 solar units, (B — V) color indices of these 

objects also span a wide range of values (O'Neil et al., 2000 |2J). All these properties indicate 
that LSB galaxies are not homogeneous by their origin and the character of evolution. 

In general, disks of LSB galaxies resemble disks of high surface brightness galaxies (HSB) 
with a scaled down density of stellar population and, usually, with lower surface density 
of gas. They have a similar character of surface brightness radial distribution (brightness 
radial profiles are close to exponential, and can also have a sharp cut-off, as for many HSB- 
galaxies). Spiral arms and bars are also observed in some of LSB galaxies (though they occur 
not so frequently as in HSB ones). The amount of young stars with high UV-luminosity per 
unit of HI mass in these systems is by an order of magnitude lower than in normal galaxies 
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(Boissier et al., 2008 [3]). 

The most essential property of LSB-galaxies is a high total mass-to-light ratio M/L 
(within several radial scalelenghts), indicating the great amount of dark matter or low mas- 
sive stars. However a high total M/L doesn't provide a clue to the location of dark mass, 
whether it is located in halo, in disk, or in both components. For clarifying the nature of 
LSB-galaxies the key role plays the problem of their disk mass. If normal stellar population 
and gas make up the bulk of a disk mass, the LSB disks have a very low surface density, and 
the total mass of these galaxies inside of optical radius must be several times higher than 
that of the baryon matter. The conclusion about the dominating dark matter contribution 
within optical radii of these galaxies is based on these arguments (see e.g. Bothun et al., 
1997 ^ and de Blok, McGaugh, 1997 [5], de Blok, Bosma, 2002 (6j, Kuzio de Naray et al, 
2008 [7J). This assumption allows one to derive the dark halo density distribution directly 
from the rotation curve of LSB-galaxies. A comparison of the obtained density profile with 
that resulted from numerical simulation of the formation of galaxies may be used to test 
different scenarios of formation of galaxies in dark halo gravitational field (see e.g. Kuzio de 
Naray et al, 2008 [7j). 

The low thickness of edge-on LSB disks, with respect to the radial scalelenght h is also 
an important argument in favor of dominating of dark halo mass contribution within several 
h and a marginal stability of their disks (Zasov, Bizyaev, 2003 |Sj, Bizyaev, Mitronova, 2009 

my 

Nevertheless, the question whether dark halo mass dominates in LSB-galaxies and how 
often it occurs remains open. For LSB-galaxies with measured rotation curves the domi- 
nating contribution of dark halo results primarily from the decomposition of the rotation 
curves, using the surface brightness and color indices estimates, which provide the infor- 
mation about the mass of stellar population. In this case it is usually assumed that the 
relation between color and M/L ratio is similar to that of HSB-galaxies, what is generally 
not so evident: disks may be significantly denser, and it will not contradict to the observed 
kinematics. 

The upper limit of a disk density may be found, using the so called maximal disk model, 
under the condition that a dark halo of a galaxy doesn't posses an inner hollow. M/Lr 
ratios of stellar population, obtained by this model for LSB-galaxies, appear to span a wide 
range of values from 1-2 to 10-16 solar units (de Blok, McGaugh, 1997 [S] , de Blok, Bosma, 
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2002 [6j, Kuzio de Naray et al, 2008 [7]), while for normal stellar population it lies in the 
interval 0,5-3 solar units. Unlike LSB-galaxies, for giant HSB-galaxies with high luminosities 
the maximal disk model gives the values of M/L, differing from those referred to the stellar 
population synthesis model usually no more than by O.ldex, what is within the uncertainties 
of photometrical method (see Kassin, de Jong, Weiner, 2006 [H]). Hence it appears that 
either a mass of LSB disks is significantly lower in comparison with that, expected for the 
maximal disk model being negligible with respect to the dark halo mass, or the disks are 
really as massive as in the HSB-galaxies and contain a large amount of invisible matter: 
low massive stars or other kinds of baryonic (e.g. cold molecular hydrogen) or non-baryonic 
dark mass, (see de Blok McGaugh 1997 j5]). 

However, the estimate of maximal disk mass is not reliable, because it strongly depends on 
the accuracy of determination of the rising part of rotation curve. A perfect examples of such 
case can be Malin-1 and NGC7589, where two different methods of the HI kinematic data 
reduction have led to the dramatically different shapes of the inner parts of rotation curves 
(F.Lelli, R.Sancisi, F.Fraternali, 2010 [12J). The disk component of rotation curve of LSB- 
galaxies usually has a peak at such distance R from the center, where the observed rotation 
curve doesn't reach its flat part yet. That's why even a small discrepancy between the radial 
scalelenghts of surface brightness and density leads to significant errors of estimation of the 
disk contribution to the rotation curve. It is worth noting, that the obtained value of a disk 
mass depends also on the assumption of the constancy of stellar population mass-to-light 
ratio. 

The alternative method of the disk mass estimation is based on the analysis of the inner 
motions or the condition of spiral structure formation in disks. This approach leads to the 
conclusion that the disks may be more massive than it follows from photometrical data. 
Masset, Bureau, 2003 [l3j have considered the numerical dynamical model of a galaxy, 
reproducing the formation of two-armed spiral structure in the extensive gaseous disk of the 
dwarf galaxy NGC 2915. Authors came to the conclusion that its disk contained dark matter, 
which distribution was proportional to that of HI mass, and was by an order of magnitude 
more massive than the latter one. Fuchs, 2003 [H] on the basis of the density wave linear 
theory has obtained for the disks of several LSB-galaxies with two-armed spiral structure 
the values of M/Lr, lying in the interval 4-16, which were significantly higher than it was 
predicted by stellar population models, but surprisingly close to the estimates of maximal 
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disk model. The existence of periodical large-scale perturbations in the disk of edge-on 
LSB-galaxy IC2233, also indicates indirectly that the disk is self-graviting (Matthews, Uson 
2008 [15]). 

In the work of Lee et al., 2004 [16] a high mass-to-light ratios of the disks of LSB-galaxies 
are explained in terms of abnormally steep initial mass function with a power law exponent of 
a ~ 3, 85 (compared with the standard Salpeter IMF a ~ 2, 35). A lower content of evolved 
stars in this case explains the lower metallicity, which, parallel with the lower amount of 
dust, leads to the lower color index of stellar population. Note that the problem with the 
interpretation of the observed chemical composition of LSB-galaxies with high values of a 
still remains (Mattsson et al, 2008 [IT]). 

The independent estimate of the upper limit of disk mass can be obtained on the basis of 
measured stellar velocity dispersion. This method, firstly applied by Bottema, 1993, 1997 
[18], [19] to HSB-galaxies, relies on the assumption of the stability of a disk against radial 
density perturbations. The denser is the disk, the higher should be the velocity dispersion of 
stars (or gas, if it makes up the bulk of disk mass) for the disk to be stable. The evaluations 
of disk mass and mass-to-light ratio for normal spiral galaxies (besides some early type 
spirals), obtained in such a way, are in good agreement with those referred to the stellar 
population synthesis models. This enables to conclude, that the disk density of galaxies is 
usually close to that, expected for the marginal gravitational stability (Zasov et al., 2004, 
2010 [20], [ZQ) 

II. THE USED METHOD 

The constraints on the disk surface density distribution may be found on the basis of 
radial velocity dispersion, using the condition of gravitational stability of stellar disk. The 
critical value of radial velocity dispersion which makes a thin, collisionless, isothermal disk 
stable against the density perturbations at all wavelengths, is usually written as: 

c cr = Qc- c T (1) 

where Q c is a Toomre stability parameter, Toomre critical radial velocity 

dispersion, cr* is a disk surface density, k is the epicyclical frequency: 

k = 2VI- y/l + (R/2Sl)(dn/dR) (2) 
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O is the angular velocity. Radial velocity dispersion can be found from the line-of-sight 
velocity dispersion measured along the slit passing through the center of a galaxy using a 
simple expression: 

c obs( r ) — ( c z • cos 2 {i) + • sin 2 {i) ■ cos 2 (a) + c r sin 2 {i) ■ sin 2 (a)) 0,5 (3) 

and the additional conditions: c r = 2Q ■ c^/k, and where c obs (r) is the line-of- 

sight velocity dispersion, c z , c^, c r are its vertical, azimuthal and radial components, i is 
the inclination of the disk of galaxy, a is the de-projected angle between the slit and the 
major axis. The coefficient k was taken to be equal to 0,6. In the model of a thin disk 
and axissymmetric perturbations, the Toomre parameter Q c — 1. But in general, the value 
of Q C (R) depends on geometrical and kinematical parameters of a disk and on the initial 
conditions of its dynamical evolution. In the current work the Toomre parameter was taken 
to vary with radius as it resulted from the dynamical modeling of 3D marginally stable disks 
(Khoperskov et al., 2003 [22]): 

Q c (r/h) = A + Ax ■ (r/h) + A 2 ■ (r/h) 2 (4) 

where Aq = 1,25, A\ = —0,19, A 2 = 0,134, and r/h - is the distance from the center of 
galaxies in the units of a disk scalelength. 

III. MASS-TO-LIGHT RATIOS 

The gravitational stability criterion can be used to calculate the upper limit of the disk 
local density at a given galactocentric distance R, and to determine the maximal contribution 
of the locally stable disk to the rotation curve. It was done for four LSB-galaxies, for which 
the distribution of the velocity dispersion and rotation curves of ionized gas and stars, 
parallel with the photometrical profiles in R-band, were obtained by Pizzella et al., 2008 [Tj 
Some parameters of the selected galaxies are listed in Table [TJ 

(1) - Name 

(2) - Distance D in Mpc 

(3) - Inclination 

(4) - Morphological type 

(5) - Absolute stellar magnitude in B-band 



(6) - Mass of HI, normalized by luminosity Lb, (parameter hie in Hyperleda [23] ) 

(7) - Notes 



TABLE I: The main properties. 



Name 


D, Mpc 


i 


Type 


M B 


hie 


Notes 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


eso-lv 1860550 


60,1 


63 


Sab 


-19,9 


1,25 




eso-lv 2060140 


60,5 


39 


SABc 


-19,65 




in group 


eso-lv 2340130 


60,9 


69 


Sbc 


-20 


1,11 


in group, in cluster 


eso-lv 4000370 


37,5 


50 


SBc 


-19,1 


1.39 


in cluster 



The criterion Q doesn't take into account a contribution of cold gaseous component to 
the condition of stability of a disk, which leads to the increasing of parameter Q c . However, 
the surface density of gas, estimated for three of four considered galaxies on the basis of 



Hyperleda [23] data (for |eso-lv2060140 such data are not available) is significantly lower than 
that of stellar disk. Taking into account the gas contribution for such ratio of stellar and gas 
surface densities following the models of Rafikov, 2001 [21] and Morozov, Khoperskov, 2005 
[25] leads to the conclusion that the overestimation of disk surface density for marginally 
stable disk does not exceed 25% of its value. Hence the influence of gaseous component can 
be neglected. 

The ionized gas rotation curves of considered LSB galaxies, obtained by Pizzella et al, 
2008 (Ij, are highly asymmetric and irregular in comparison with the stellar ones, in this 
connection these authors recommend to use the stellar rotation curves in dynamical mod- 
eling. But stellar rotation curves are not as extensive as gas ones, so both stellar and gas 
rotation curves were used in the current work to construct the dynamical models. Stellar 
rotation curves were corrected for asymmetric drift according to Binney, Tremaine, 1987 
[26]: 

v \ = v 2 r + c 2 . (o, 5 • d{lnv r )/d(lnr) -0,5 + r/h- d(lnc 2 ,) / 'd{lnr)) , (5) 

where v c and v r are circular and observed rotation velocities correspondingly, c r is the radial 
velocity dispersion, h is radial scalelenght of the disk. 
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Surface density radial profiles of marginally stable disks, obtained by applying the local 
criterion (JTJ, were used to estimate (M/L R ) disk for different galactocentric distances (see 
Fig.l). The filled and open symbols in Fig.l correspond to the gas and stellar circular 
velocities. The dashed lines denote (M / Ln)di S k ratios, predicted by stellar population models 
(Bell, de Jong, 2001 [27]) for the metallicity Z = 0,008, the solid lines correspond to the 
mean values of (M/L R )di S k, which were accepted for the decomposition of rotation curves. 



ESO-LV 1860550 



R, kpc 



ESO-LV 2060140 




R, kpc 



ESO-LV 2340130 



R, kpc 



ESO-LV 4000370 



R, kpc 



FIG. 1: Radial profiles of {M / L R )di s k ratios, estimated from the local gravitational stability crite- 
rion. Open and filled circles correspond to the gas and stellar circular velocities. The dashed lines 
denote {M/L R )di s k ratios, predicted by stellar population models (Bell, de Jong, 2001 [27]). The 
solid lines correspond to the mean values of (M / L^jdiski which were used for the decomposition 
of rotation curves. 

From Fig. 1 it follows that in spite of some differences, the distributions of (M/Ln)di S k, 
obtained on the basis of both gas and stellar rotation curves in three of four cases are nearly 



constant (for eso-lv2060140 this ratio is roughly constant only in the inner part R < 5 
kpc). Central parts of radial profiles of (M/L R ) d i S k of the galaxies with significant bulges 
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were excluded from consideration. From Fig. 1 it can be seen, that (M / L^jdisk ratios, 
corresponding to the stellar disks marginal stability criterion, are several times higher than 
those predicted by stellar population synthesis models. It is worth reminding, however, 
that the marginal stability criterion allows one to obtain only the upper limit of disk mass. 
Nevertheless, the approximate constancy of (M / L R )di S k{R) for marginally stable disks gives 
evidence that this model may be acceptable for real galaxies. 

IV. THE ROTATION CURVES MODELING 

The circular rotation curves of gas and stars were decomposed into the separate com- 
ponents: the exponential disk, the pseudo-isothermal dark halo and, if necessary, King's 
bulge using the values of (M/Lr)^^, found from the gravitational stability condition (see 
Fig. 1 and Table |ll|). The radial scalelengths of disks were taken to be equal to the photo- 
metrical scalelengths, obtained in R-band by Pizzella et al., 2008 pp. The decomposition 
of the rotation curves of stars after correction for the asymmetric drift (filled circles) and 
gas (open circles) is illustrated in Fig. 2. For |eso-lv 2060140 two modeled rotation curves, 
corresponding to two different [Mj L R ) disk are presented, since one can't unambiguously 
determine (M/L R ) disk for this galaxy. 

From Fig. 2 it follows that the models, corresponding to the marginal gravitational 
stability criterion, provide a reasonable fits to the observed rotation curves. In Table |TT] the 
following results of rotation curve modeling are listed (the masses of the components within 
R = 4h were found by means of the extrapolation of the rotation curves to this radius): 

(1) - Name 

(2) - (M/L R ) disk ratio, which was used for the modeling 

(3) - The disk mass inside R = 4h 

(4) - The total mass inside R = Ah 

(5) - The dark halo mass inside R = 4h 

As it follows from Table [IT] the obtained mass of the halo for all considered galaxies is 
no more than twice as high as that of the disk within R = Ah. These results indicate that 
the disks of LSB-galaxies considered here might be significantly more massive, than it can 
be expected from their low brightnesses for normal stellar population. The disk masses of 
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R, kpc 




FIG. 2: The decomposition of the rotation curves of ionized gas (filled circles) and stars (open 
circles) 

different LSB-galaxies, obtained in this work (filled circles) and by Fuchs, 2003 [H] (open 
circles) are compared with those referred to the maximal disk models in Fig.3. From Fig. 3 
it follows that the mass models, obtained in the current work and by Fuchs, 2003 [14J are 
close to the models of maximal disk, where stellar disk has maximal possible contribution 
to the rotation curve. 



10 



TABLE II: The results of modeling 



Name 
(1) 


(M/L R ) disk 
(2) 


M disk • W 10 M sun 
(3) 


M tot ■ W w M sun 
(4) 


M hal0 ■ W w M sun 
(5) 


eso-lv 1860550 


7,5 


6,2 


13 


6,9 


eso-lv 2060140 


8,0 


7,0 


12 


4,9 


eso-lv 2060140 


8,8 


7,7 


8,4 


0,7 


eso-lv 2340130 


5,7 


4,7 


9,4 


4,1 


eso-lv 4000370 


9 


5,0 


7,5 


2,3 




FIG. 3: Comparison of LSB disk masses, obtained using my models (filled circles) and models 
of Fuchs, 2003 p3] (open circles) with those, referred to the maximal disk models. Black line 
corresponds to M disk = (M disk ) 

max' 

As it was mentioned above, LSB-galaxies are often used to test the models of galaxy 
formation in the gravitational field of a dark halo by comparing the observed density profiles 
of these objects with those predicted in terms of cosmological modeling, neglecting the 
contribution of disk component. The results described above indicate that this approach 
may not work at least for considered LSB-galaxies. Moreover, according to Kuzio de Naray 
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et al., 2008 [7], the density profile predicted by cosmological modeling (NFW profile) badly 
enters in the maximal disk models of LSB-galaxies (in other words, the NFW and "max disk" 
models are incompatible). In this case, if the disks of the galaxies considered here aren't 
dynamically overheated, and their mass models are really close to the models of maximal 
disk, it can be an argument against the universality of NFW profile. 



V. ON THE DENSITIES AND THICKNESSES OF STELLAR DISKS 

The stellar velocity dispersion data allows to obtain not only disk surface densities (or 
their upper limits), but also half-thicknesses of stellar disks (or their lower limits), cor- 
responding to their stability condition. A half-thicknesses of isothermal self-gravitating, 

2 

marginally stable disks: 2* = at the fixed distance R = h may be compared with 

those calculated using the densities, resulted from the photometrical data. In the later case 
a self-consistent equilibrium equations were solved for stellar and gaseous components of a 
disk in the gravitational field of rigid massive spheric pseudo-isothermal dark halo. 

A comparison of the disk parameters, obtained using marginal stability condition 
(column a) with those referred to the stellar population modeling (column b) is presented 
in Table [TTT1 

(1) Name 

(2) - z*/h ratio at the galactocentric distance R = h 

(3) - Central surface density of stellar disks <7o* 

(4) The ratio of halo and disk masses inside R = 4h 

(5) - Radial scalelenghts of the disks 



For eso-lv2060140 in case (a) the value (M/L^disk — 8 wa s used. From Table III 



it 



can be seen, that the photometric models lead to the higher contribution of dark halo to 
the rotation curve in comparison with marginally stable model. The disk surface densities, 
estimated for normal stellar population models give 2-2,5 times thicker disks in comparison 
with the marginal stable ones. The estimates of z*j h, corresponding to the marginal stability 
criterion, better agree with those expected for LSB galaxies from the relationship between 
the deprojected surface brightness and disk thickness, found by Bizyaev, Mitronova, 2009 
[9], where the objects with lower surface brightness have lower z*/h. 
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TABLE III: A comparison of parameters estimated from gravitation stability criterion and pho- 
tometric models. 



Name 
(1) 


z*/h, pc 
(2) 


0"o*, M sun /pc 2 
(3) 


Mhalo/M disk 

(4) 


h, kpc 
(5) 




a 


b 


a 


b 


a 


b 




eso-lv 1860550 




0,1 


0,2 


1000 


240 


1,1 


11 


3,6 


eso-lv 2060140 




0,2 


0,4 


450 


68 


0,7 


9,1 


5,3 


eso-lv 2340130 




0,1 


0,3 


590 


150 


0,9 


7,4 


3,8 


eso-lv 4000370 




0,2 


0,4 


520 


85 


0,5 


26 


4,1 



VI. CONCLUSIONS 



The disk masses of four LSB-galaxies, obtained using the local gravitational stability 
criterion, appear to be significantly higher, than it can be expected from their photometrical 
data in the assumption of normal stellar population. However they agree with the observed 
rotation curves. The corresponding mass models are close to the models of maximal disk. 
Hence, either the disks of these galaxies are strongly dynamically overheated and thick, or 
their densities and masses are quite typical for HSB-galaxies. In this case their high mass- 
to-light ratios can indicate that either they contain a large amount of low massive stars (a 
heavy bottom IMF), or a significant fraction of dark matter is contained in the disks, either 
in baryonic (e.g. cold molecular gas, see Pfenniger, Combes, 1994 [30]) or non-baryonic 
form. 

It seems that at least some LSB-galaxies are not dark halo dominated systems which 
possess the disks of low density. In this case LSB galaxies differ from HSB ones not by the 
high contribution of dark halo, but rather by peculiarities of their disks in the presence of 
dark halo of moderate mass. 
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